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ABSTRACT 


This thesis aims at modeling of the smart actuators and dicusses their application in 
reduction of vibrations in a highly simplified helicopter model. A smart actuator is 
a stack of alternate layers of metal and piezoelectric material subjected to externally 
applied force and/or electric potential. The deflection of such an actuator is dependent 
on these externally applied forces and voltages. A finite element modeling is carried 
out in the present study l.o formulate relation between deflection and the externally 
applied forces and voltages. When there is a cyclic variation of electric field, piezoelectric 
materials exhibit polarization-electric field (P-E) and strain-electric field (e-E) hysteresis 
effects. These curves are commonly termed as “ hysteresis loo{)” (for (P-E) variation) and 
“ butterfly loop” (for f-E variation). The present formulation tries to capture both these 
effecls by following two different approaches, namely the complex stiffness model and the 
domain switching model. The domain switching model has proved to be a much preferred 
model as com|)ar('d to the otlier. Hysteresis loop is a representation of the energy loss in 
the system. This energy loss heats up the i)iezoelcctric material. A heat transfer analysis 
has been carried out to find the variation in temperature distribution within the actuator 
with change in excitation frequency of the system. Tlie effecd, of including convection mode 
of heat dissipation is emphasized througli this analysis. A highly simplified helicopter 
model is considered as a two degree of freedom system representing hub and fuselage 
connected by a S[)ring-dampcr-pi(;zo-actuator-scnsor system, and a closed loop active 
vibration control is carried out on this model. In this study, the model consisting of 
alternate layers of metal and pioezoelectric material is used as a sensor as well as an 
actuator. A comparison between the passive and active vibration isolation is highlighted. 
It was shown that the model developed minimizes the vibratory loads of both the rotor 
and fuselage significantly. 
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NOMENCLATURE 


<1, b, aj, aj, bj, C(, C2 Consfcanf.s 
d, b Vibration absorber dimensions 

A Area 

c StifFness cocfTieicnts 

C Dami)ing 

Critical damping coefficient 
( 1 , e, q, g Piezoelectric coupling coefficient 
D Electric disi)lacement 

dV Elemental volume 

dS Elemental surface 

c' Piezoelectric stress coefficient 

E Electric field 

E Young’s modulus for metal layer 

f Frequency in Hertz 

F Force 

Fq, Fy<^ Forces al. the al.tachment points (i.e., at hub-gearbox and fuselage-gear box 
attachment) 
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Fpc, Fp^ Forces applied on the sensor/ actuator 
h Heat transfer coeflicient 

k Thermal conductivity 

^M’ Thermal conductivity of metal and piezo material respectively 

K Stiffness 

Lpp Lp Length of metal and piezo layer rcsi)cctivcly 
Mpj, Mp Mass of gear box/hub and fuselage rcsjrectively 

N Effective number of aligned dipoles 

P Electric polarization 

P Perimeter 


Q 

<')Q 


R 
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AT 


Heat flux 

Heat generated/unit cycle/nnit volume 
Heat generated 
Heat flow 

Variation in (piantity of heat 
Compliance 
Curie Temperature 
Temperature (T the bar 
Ambient tem[)crature 
Temperatures at different hiyers 
Time period 
Change in temperature 
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Vy\^ Actuation voltage 
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Greek Symbols 

Of Coefficient of expansion 
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jj Dielectric susceptibility 

Lj Frequency in rad/s 

ujref Reference frecpiency 

(/) Phase angle 

?/; Angle in actuator/sensor tnechanisin 

T Electric potential 

Induced potentia.1 
a Stress 

0 Temperature difference 
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Chapter 1 


INTRODUCTION 


1.1 Overview 

"[’ho technological developments in acrosjiace engineering have created new avenues of 
research particularly in the areas of structural health monitoring, vibration control, shape 
and noise control using the concei)t of smart or intelligent structures. Under external 
loading conventional structures deform in a passive manner with no control over the 
final deformed configuration. Whereas in smart structures, tlic structural deformation is 
contiimously monitored and controlled (using the distributed sensors and actuators) to 
achieve the desired configuration of the deformed state of the structure. The concept of 
a smart structure can be well (hiscribed as "a sjfslcm which ha.’i inl.rvisic sensor, actuator 
and control mechanisms whereby it is capable of sensing a stirnuhi,s, responding to it, and 
reverting to its original state after the stimuhi.s is jxanoved". 'The stimulus may corrcs|)oad 
to a stress, strain, light, electric field etc., whereas the response represents the motion 
or change in optical properties, modulus, surface tension etc. Smart structures can be 
basically classified as: 

1. Passively smart 

Structures which have the ability to respond to a stimulus in a useful manner, 
without assistance of electronic controls or feedback systems 

2. Actively smart 

Structures which utilize a feedback loop to accelerate the recognition and response 
process 
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3. Very smart or' intelligent 

Structures which utilize the nonlinear properties of the sensor, actuator, memory 
and/or feedback systems to tune the response behavior. 

A smart strucl.urc comprises of three main components and tliey are;- 

1. Sensor 

To monitor environmental changes and generate signals proportional to the chang- 
ing measurand. Sensor applications are based on the direct piezoelectric effect 
(strain-induced electric field). 

2. Actuator 

The actuators are used to change the properties of smart structure in order to 
achieve the desired response. The actuator applications are based on the indirect 
piezoelectric effect (apj)lied external voltage-induced strain). 

3. Control systems 

The control S 3 'stem continually monitors the sensors signal, processes the informa.- 
tion in order to determine the required action . 

The distributed sensing and actuation can be implemented by embedding piezoelectric 
material in the structure. 

1.2 Piezoelectric Materials: History and Processing 

r’i(!zoele(;tricity stcuns from the Greek word piezo for i)ressur(!. It follows tliat a piezoelec- 
tric material develops a potential across its boundaries when subjected to a mechanical 
stress (or |)ressuie), and vice versa, when an electric field is applied to the material, 
a mechanical deformation ensues. Piezoelectric materials therefore fall in the class of 
smart materials. Perro-electricily is a subgroup of piezoelectricity, where a spontaneous 
polarization exists that can be reoriented by application of an alternating electric field. 

1.2.1 History 

Piezoelectricity is a pro[)er(.y possessed by a group of inaterials, discovered in 1880 fry 
Pierre and .fac(iu('s Curie during their study of i,he elfects of pressure on the generation 
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of electrical charge by crystals such as Quartz, Tournialiue, and Rochelle salt. In 1881, 
the term "piezoelectricity" was fust suggested by W.Hankel, and the converse effect was 
deduced by Lipniann from thermodynamics principles. In the next tliree decades, collab- 
orations within the European scientific community established the field of piezoelectricity; 
and by 1910, Voigt’s "Lerbuch der Kristallphysic" was published and became a standard 
reference work detailing the complex elcctro-thcrmo-elastic relationship in piezoelectric 
crystals. However, the complexity of the science of i)iezoclcctricity made it difficult for 
it to mature to an application until a few years later. Langevin et al. developed a piezo- 
electric ultrasonic transducer during World War I. Their success opened up opportunities 
for piezoelectric materials in undcrwai.er applications such as ultrasonic transducers, mi- 
crophones, accelerometers, etc.. In 1935, Busch and Scherrer discovered piezoelectricity 
in potassium dihydrogen phosphate (KDP). The KDP family was the first major family 
of piezoelectrics and ferro-electrics to be discovered. 

During World War II, research in piezoelectric materials expanded to the U.S., the 
Soviet Union and .Ta,pan. Up until then, limited performance by these materials inhibited 
commercialization but that changed when a major breakthrough came with the discovery 
of barium titanato and lead zirconate titanatc (PZT) in the 1940s and 1950s respectively. 
These families of juatcrials exhibited very high dicicci.ric and piezoelectric properties. 
Furthermore, they offered the {)ossibility of tailoring their behavior to specific responses 
and apj)lications by the use of doj)anf,’s. To date, PZl’ is one of the most widely used 
pi('Z()eIectric materials. It is noted that most commercially available ceramics (such as 
Ijarium titanatc and PZT) are based on the perovskitc structure (Fig. 1.1). The per- 
ovskite structure (ABO 3 ) is' the simi)lest arrangement where the corner-sharing oxygen 
octahedra are linked together in a regtdar crd)ic array with smaller cations (Ti, Zr, Sn, Nb 
etc.) occiit)ying the central octahedral B-sitc, and larger cations (Pb, Ba, Sr, Ca, Na etc.) 
filling the interstices between octahedra in the larger A-sitc. Compounds such as BaTi 03 , 
PbTi 03 ,PbZr 03 , NaNbOg and KNbO;} have Ireen studied at length and their high tem- 
perature ferroelectric and anti ferroelectric phases have Ijeeu extensively exploited. This 
structure also allows for multiple substitutions on the A-site and B-site resulting in a num- 
ber of useful though more com])lex compounds such as (Ba,Sr)Ti 03 ,(Pb,Sr)(Zr,Ti) 03 , 
Pb(Fe,Ta) 03 , (KBi)Ti 03 etc. 

Starting around 19G5, several .Japanese companies focused on developing new pro- 
ces.ses and ai)plications, and opening new commercial markets for ]>iezoclcctric devices. 
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The success of flic Japanese effort attracted other nations, and today the needs and uses 
extend from medical applications to the communications field to military applications 
and the automotive field. 

A review of the early liistory of piezoelectricity is found in the work of W.G.Cady 
(Ref. [11). 

1.2.2 Piezoelectric Ceramic Processing 

The fabrication of most bulk piezoelectric ceramic starts with powder preparation. The 
powder is then pressed to the required shapes and sizes, and the grain shapes are in turn 
processed to mechanically strong and dense ceramics. The more important processes 
that influence the iiroduct characteristics and properties are ])Owder preparation, powder 
calcining and sintering. The next steps are machining, electroding and poling: application 
of a DC field to orient the di|)oIes and induce piezoelectricity. 

The most common powder {ueparation is the mixed oxide route. In this process, pow- 
der is prepared form the appropriate stoichiometric mixture of the constituents oxides. 
In the case of lead zirconatc titanate (PZT): lead oxide , titanium oxide, and zirconium 
oxide are the main compounds. Depending on application, various dopant’s are used to 
tailor the properties of interest. PZT ceramics arc rarely utilized without the addition 
of dopant’s to modify some of their properties. A-sitc additives tend to lower the dissi- 
pation factor, which affects hea,t generation, bid. also lowers the piezoelectric coefficients; 
for this reason I, hey are mosl.ly used in ultrasonics and other high frexiuency a.pplications. 
B-site dopants increase the iiiezoelectric coefficients but also increase the dielectric con- 
stant and loss. They arc utilized as actuators in vibration and noise control, benders, 
optical positioning application etc. Figure 1.2 shows a flowchart of the mixed oxide route 
for making PZT ceramics. Mixing of the powders can be done by dry-ball or wet-ball 
milling, both methods having advantages and disadvantages: wet ball-milling is faster 
than dry-ball milling; however, the disadvantage is the added step of liquid removal. The 
most common method for making PZT ceramics is through wet-ball milling; ethanol and 
stabilized zirconia media are added for a wet milling process. A vibratory mill may be 
used rather than a conventional ball mill; it w^is shown by Ilerner (Ref. [2]) that this 
process reduces the risk of contamination by the balls and the jar. Zirconia media are 
used to further reduce the conta.mination risks. The calcination stei) is a very crucial step 
in the processing of PZT ceramics; it is important that the crystallization be complete 
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and that the perovskitc pliasc forms during this si.ep. The goals are to remove any or- 
ganics, water or other volatiles left after mixing; to react the oxides to form the desired 
phase composition before the ceramic is j)roccssed to useful devices; and to reduce volume 
shrinkage and allow for better homogeneity during and after sintering. 

After calcining, a binder is added to the powder, then the mixture is shaped usually 
by dry-pressing in a die for simple shapes, or extrusion, or casting for more compli- 
cated bodies. Next, the shapes are sintered: placed in an oven for binder burn-out and 
densification. 

The major problem in the sintering of the PZT ceramic is the volatility of PbO at 
about 800® C. do minimize this problem, the PZT samples are sintered in the presence 
of a lead source, such as PbZrOa, and placed in closed crucibles. The saturation of the 
sintering atmosphere with PbO minimizes lead loss from the PZT bodies. Sintering can 
now be carried out at temperatures varying between 1200-1300® C. Despite precautions, 
there is usually a resulting loss of 2%-3% of the initial lead content. 

After cutting and machining into desired sha|)es, electrodes are applied and a strong 
DC field is used to orient the domains in the polycrystalline ceramic. DC poling can 
1)0 done at room temperature or at higher temperatures dej)ending on the material and 
the composition. Tlie i)oling process only partially aligns the dipoles in a polycrystalline 
ceramic, and the resulting polarization is lower than that for single crystals. 

Tiiis processing tcchni(iue [uesent many uncerl.aititles and the presence of a wide 
number of other faljrication techniques is an indication that there is a great need for the 
production of reliable PZT ceramics with optimum properties and micro structure. One 
problem often encountered is the deviation from stoichiometry. This problem is often due 
to impurities present in the raw materials as well as the lead loss during the sintering 
process, which invariably residts in substantial alterations of the PZT properties. As a 
result, the clastic properties can vary as much as 5%, the piezoelectric properties 10% 
and the dielectric properties 20% within the same batch (R.ef. [2]). Also, the piezoelectric 
and dielectric properties generally suffer if there is any lack of homogeneity due to poor , 
mixing. It is imirortant then that the constituent oxides be intimately mixed. In the 
metliod descrilred above, however, the constituents are solid solutions and it has been 
shown that an intimate mixing of solid solutions is difficult if not impossible. More 
information on the preparation of [)iezoclectnc ceramics can be found in Jaffc, Cook and 
•Taffe (Ref. [3]). 
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1.2.3 Definition 


Piezoelectricity is a physical phenoTncnon which is related to the microscopic structure of 
the solid. Some ceramic materials become electrically polarized when they arc strained; 
this phenomenon is rcfci red to as direct piezoele.clr ir. effect. The direct piezoelectric elTect 
is always accompanied by tlie converse piezoelectric effect where a solid becomes strained 
when placed in an electric field. The microscopic origin of the piezoelectric effect is the 
displacement of ionic charges within a crystal structure. In the absence of external strain, 
the charge distribution williin the crystal is symmetric and the net electric dipole moment 
is zero. However when an external stress is applied, the charges are displaced and the 
charge distribution is no longer symmetric. A net polarization develops and results in an 
internal electric field. A material can only be piezoelectric if the unit cell has no center 
of inversion. 

The interaclion among the three fundamental constituents namely mechanical, elec- 
trical and thermal loading on piezoelectric ma.teria.ls is shown in Fig. 1.3. The various 
interactions that basically exist can be classified as thermo-elasticity, piezoelectricity and 
pyroelectricity (lief. [Ij and (4|). An electric field E causes j)iezoelectric stress a = e'E, 
where e' is the a[)pro[)riate piezoelectric stress coefficient. Similarly a strain e causes the 
electric polarization F = cc. In a similar manner, the pyroelectric constant p relates a 
change in temperature AT with polarization P; the arrow from E to 5Q indicates the 
electro-caloric effect (i.e., the variation in the quantity of heat Q on application of an 
electric field E), which is usually expressed as a relation between the change in temper- 
ature Al" and E. The l.hcrmal coefficient of expansion a; relates AT with e; the arrow 
from (7 to 6Q with a coefficient b relating them indicates the therrnoclastic effect. Stress 
a and strain c arc', related by the ecpiation of the form e- = sa. In an analogous manner, 
polarization P and electric field E are related by dielectric susceptibility r/; and quantity 
of heat 5Q is related to change in temperature AT through specific heat c. The arrows 
indicate the directions in which the various effects usually take place. The additional 
arrow between P and E indicates that an electric field may exisl. l>y virtue of polarization . 
in the material. 

Each of the nine straight linos in Fig. 1.3 usually represent the primary or “true effeefP. 
However, in every case there is at least one other jjath over which the process can take 
place, unless the coefficieni, are zero for I, hat particular process for the material (or crystal) 
under consirleration. The roundabout effects are usually denoted as “false or secondary 
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effects”. An example is the false pyroelectric effect cine to piezoelectric action, which may 
be several times greater than the direct effect itself. The i)riinary pyroelectric effect is 
indicated by the path AT — > P, while the secondary effect follows the path AT — > e — )• T. 
The reason for this complication is dnc to the fact that every pyroelectric crystal is also 
piezoelectric. Therefore a cliangc in temperature of a,n unconstrained crystal causes a 
deformation and this in turn produces the secondary i)olarization of piezoelectric origin 
superposed on the primary pyroelectric polarization. Therefore it is evident that primary 
pyroelectricity would be observed only in a conrjdetely clamped crystal. 

In a similar maimer when an elastic compliance coefficient s is measured by observa- 
tions of a and e a piezoelectric polarization P is produced, which if the crystal is not short 
circuited gives rise to an electric field E, which in turn modifies the value of stress a. In 
addition, through thermoelastic effect there is a change in the temperature which affects 
the value of strain e. Hence it becomes important to specif}'’ the conditions under which 
the material coefficients arc measured and they must be properly incorporated in the 
analysis of smart structures^ having a base structure with embedded piezo layers. Under 
cyclic variation of electric field, piezoelectric material exhibits a hysteretic (polarization- 
electric field (P-E)) effect leading to energy loss. This energy goes in for developing heat 
i.e., hysteresis loss corresponds to heat develojred/dissipated. As a result, the actuator 
heats up tremendously and this sometimes becomes the cause for the damage of the actu- 
ator, In addition to the P-E hysteresis the polarization switcliing leads to a strain-electric 
field hysteresis termed as “buttcrily loop”, fioth these loops are shown in Figs. 1.4, and 
1 . 5 . 

It is due to the direct and converse effects l.hal, the jiiezoclectric material embedded 
smart sf.ructurcs find great use as both sensors and acl.ual.ors. Smart actuators having 
piezo stack have Ireen developed for generating active control forces in active vibration 
isolation systems. 

1.3 Literature survey 

The literature survey has been carried out in three distinct directions of research, namely: 

• One dealing with the modeling of the smart actuator. 

• Formulation of a model which would general.e hy.sl.crc,sis iind hutierfly loops. 
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• Vihnitioii control in liclicopters using smart actuator. 

1.3.1 Modeling of Smart Actuator 

The authors of ftof. |2] have provided a detailed review of luezoclcctric ceramics analysis 
and characterizatioti. The foctis is on polycrystalline ceramics. To thoroughly grasp the 
l)ehavior of a piezoelectric [)olycrystaIline ceramic, a basic understanding of the ceramic 
itself should not be overlooked. I'his review presenl.s a brief introduction of the history 
of piezoelectricity and a discussion on j)roccssing of the ceramic and development of the 
constitutive relations that define the behavior of a piezoelectric material. It also covers 
the most common methods of measuring the parameters of interest in piezoelectricity. 

In Ref. |5| Venkatesan and Upadhyay presented a conpuchensive description of the 
jihysical interactions in a piezoelectric material used in development of smart structures. 
The governing electro-thermo-elastic eejuations for the analysis of smart structures have 
been developed from the fundamental laws of conservation and thermodynamic principles. 
Further in Refs. [4], [G] and [7|, using a finite element model, the deformation of a smart 
beam has been analyzed for actuation and sensing cases, by solving the coupled electro- 
clastic ijroblem. The key dilfercnce lietween sensing and .actuation has been brought out 
by solving for the electric potential and the di,sj>laccmcnt vari.ables in the piezo materials. 
It is also brought out that the interaction between polarization and electric field in the 
piezo continuiuu Ic.ads to tion-linear distributed body forc'c .and non-symmetric stress 
tensor. These non-linear cficcts arc found to have significant influence on the deformation 
of a smart beam under actuation. 

Sirohi and Choijr.a (Ref. [S|) investigated the behavior of piezoccr.amic actuators 
under different l.ypes of excil.ations and mechanical lo,a,ding. The research focuses on the 
aijplication of these .actuators to the develoinnent of smart rotor systems. The free strain 
response of the actuators under DC excitation is investigated experimentally along with 
the associated drift of the strain over time. Effect of tensile stresses on the DC respon.se is 
investigated. Strain response to an AC excitation is also investigated for a free actuator 
and for a pair of .actuators surface bonded to a host structure. The phase lag of the strain 
response and non-linear hysteretic effects have Ireen oljserved and the power consumption 
h.as been validated by the impedance method. Additionally, depoling of the .actuators is 
discu.ssed, along with the fe.a.sibility of recovering perform.ance by repoling ip the event of 
accidental depoling. 
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1.3.2 Modeling of Hysteresis and Butterfly Loops 


A phenomenological constitutive law for feiroelectric switching due to multi-axial me- 
chanical and electrical loading of a polycrystalline material was developed by Robert 
and Chad (Ref. [9]). The constitutive law succe-ssfully reproduces hysteresis and but- 
terfly loops for ferroelectric ceramics. The pai)er provides a relatively simple law that 
can simulate the nonlinear switching behavior of polycrystalline ferroelectrics subject to 
multi-axial stress and an arbitrary electric field. 

In Refs. [lOj and [11| an attempt is made to develop a mathematical model for the 
hysteresis curve by employing a hysteretic restoring force model (Refs. [12] and [13]) 
to capture the effects due to change in temperature, pressure and electric field. The 
parameters representing the model have been idejitiliod and their influences on the nature 
of hysteresis curve ha.ve been analyzed. It has been seen that there are strong interactions 
among temperature, pressure and amplitude of electric field. A small change in one factor 
can alter the behavior due to other effects. 

Chen and Monl.gomery (R.cf. |14j) present a macroscopic theory leading to the hys- 
teresis and butterfly loop in ferroelectricity. This theory is based on the notion that as 
domains switch under the action of an external electric field the number of dipoles aligned 
in the direction of the field also alters and it obeys a rate law which depends non-Iinearly 
on the field. In addition, the paper specifies the simplest constitutive relations of the 
stress and electric displacement which together with the rate law are sufficient to yield 
the hysteresis looi) and the butterfly loop. 

A phenomenological forimdation of polarization reversal of piezoelectric materials is 
proposed based on the dynamics of domain switching in Ref. [15] . This formulation pro- 
vides a method to describe the hysteresis in i)iezoelecti icity as well as in electro-magnetics. 
It is shown that, a good approach to ilescribe the nonlinear indiiced strain field behav- 
ior and electromechanical hysteresis in piezoelectricity is by combining the macroscopic 
phenomenological asi)ects with the microscopic material proj^erties. A one dimensional 
f,hermo-electro-mechanicaI constitutive model for piezo ceramics which undergo polariza- 
tion reversal is {)rcscnted using a continuum mechanics approach. This model is based 
on thermodynamic j)rinciples and reflects the essence of the electromechanical behavior 
of piezoceramics siniply and reasonably well. 

A model that bridges micro- and macroscopic couccijts was developed by Narayan, Lu 
and Ilanagud (Ref. (IGj) to describe the nonlinear material and thermodynamic bchav- 
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iols of piezoclocfiic cernniics. The doinaiii evolufiotis was described by the evolution of 
tlic associated internal varialdcs, which are the parameters of the assumed domain orien- 
tation distribution. The bulk material behavior w'as obtained by simply using the domain 
distribution to collect all the domain behaviors. Using this model, the thermodynamic 
parameters wme evaluated. 

1.3.3 Vibration Control using Smart Actuators 

The problem of the formulation and solution of an active vibration control scheme for 
helicopters, based on the concept of active control of structural response is addressed 
in Ref. [17j. First, using a mathematical procedure employing the Fischer information 
matrix, optimum sensor locations have been identified in a three-dimensional model of 
a flexible fuselage structure. Then, using the measurements from these optimal sensor 
locations, a multi-iii])ut/multi-output control probhnn has been formulated and solved 
to obtain the acl.ive control foices rcKiuired for vibration minimization in the helicopter 
fuselage. It is observed that sen.sor locations have a significant influence on the level 
of vibration reduction in a fuselage structure. A leview on the sensor placement in 
distributed parameter systems and their influence on vibration reduction can be found 
in Refs. [18|-{21j . 

The equations of motion for the coupled rotor/isolator/fuselage dynamical system 
are formulated in Ref. (22| for the prediction of viluation in helicopters during forward 
flight. Using a combined Flociuet theory/fre(iucncy-responsc technique, the equations 
were solved to predict the vibratory loads at the various locations of the helicopter. In 
addition, the influence of the isolator on the hub and blade loads was established. It 
was shown that the inclusion of an isolator between the hub and the fuselage reduces the 
fuselage vibral.ory loads and increases the hub loads significantly. 

The problem of vibraUon reduction in helicopter fuselages using the active control of 
structural response (ACSR) is addrc.ssed in Ref. [231. First a balance-realization-based 
order reduction is employed f.o reduce the size of the problem. It was observed that 
contrary to open-loop control, the technique of closed-loop control reduces the vibratory 
levels both in the fuselage and the gearbox. 
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1.4 Objectives 

The main object, ivcs of (,hc present, study are ; 

1. Modeling of hysteresis o,nd butterfly loops : To forumlate a iriodel that would 
generate pohirization-electric field {hysteresis) and strain-electric field (butterfly) 
loops 

2. Modeling of smart actuator : To analyze the deformation and induced potential 
developed in a smart aet.uator using a finite element model, by solving the coupled 
electro-elastic problem of a stack actuator. 

3. Heat transfer analysis of a smart bar : To study the heat transfer in a smart 
actuator taking into consideration that heat dissipation to the surroundings happens 
through the plienomena of conduction and convection. The importance of including 
convection is highlighted. 

4. Vibration control using a active vibration isolation concept employing posi- 
tion feedback control. 

1.5 Organization of tlie Thesis 

In Chapter 2, the constitutive equatio?is requirc<l to generate the hysteresis and butterfly 
loops are generatcfl. Two approaches are used namely: linear complex stiffness approach 
and the other is the improved model which incorporates (,hc concept of domain switching. 
The relative difference in the two models at various operating frcriucncies and electric 
fields is analyzed. 

Chapter 3 introduces the study of heat transfer in the smart actuator. The basic heat 
traiisfer phenomena i.e., conduction and convection are considered while formulating the 
erjuations. The importance of convection heat transfer in the formulation is brought out. 
The study is carried out at various frequencies and electric fields to find the temperature 
distribution in smart actuators, wliich is useful for practical real-life applications. 

Chapter 4 deals with the dcscrii)tion and formulation of a layer-by-layer finite ele- 
ment analysis for smart stack actuators. The problems of sensing and actuation have 
been addi'cssed here. The distribution of induced potential in the piezo material duo to 
externally applied potential and due to mechanical straining has been studied. 
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In Chapter 5 the problem of vibration control us addressed. A highly simplified model 
of a helicopter is considered fxs a two degree of freedom system representing hub and 
fuselage connected by a spring-damper-piezo-acf.uator-sensor system. It is shown that 
the inclusion of active vil)ration isolator between the hub a.nd the fuselage reduces not 
only the fuselage vibratory loads as well as the the hid) loads significantly, contrary to 
the passive isolation systein. 
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MIXING of OXIDES: 



Figtirc 1.2; Mixcd-Oxiclo-Ronte of i)rni)ariiig PZT 
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Figure 1.5: Butterfly loop 


Chapter 2 


MODELING OF HYSTERESIS AND 
BUTTERFLY LOOPS 


A fcrrooloctric inatei ial possesses spojil.aneous po!a.i i7,a.ti()ii tlia,t can be reversed in direc- 
tion by application of a realizable electric field. This phenomenon is observed over some 
temperature range. Most ferroelectric materials have a Curie temperature, T(j below 
which they are irolar and above which they are not. Ap[)lyiiig a large alternating electric 
field causes the j)olarization to reverse, and this gives rise to the ferroelectric hysteresis 
loop, relating the i)olarization P to the electric field E, The relation between the electric 
displacement D, polarization P and electric field is given by the linear equation 

(2.1) 

I o 

where So is the j)crmittivity of free space (--8.85*10'^'^ C/V-rti). A typical hysteresis 
loop is illustrated in Fig. 2.1. Let us begin with a crystal with randomly oriented dipoles 
so that there is no net polarization (i.e., average polarization is zero). When an electric 
field is applied, the dipoles begin to line up with the field points 1-2. As the field is 
increased, eventually all the dipoles are aligned in the direction of the field and maximum 
polarization (also denoted as saturation jrolarization) Pg^q, is obtained (point 2). When 
the electric field is removed, a remanent polarization denoted by point 3 remains due to 
interaction between dipoles and therefore the material is permanently polarized. When 
the direction of electric field is changed, the dipoles also change their direction. The 
coc’.rcive field h]c (point 4) removes the polarization and randomizes the dipole directions. 
If the field is inenursed further, saturation occurs vvith opposite polarization. As the field 
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coiitimics to alternate, a hysteresis loop is formed and the area of the loop is related to 
the energy required to switch the polarization from one direction to the other. 

The hysterctic liehavior of the piezoceramics is a complex issue, sensitive to applied 
stie.s.scs and temperature. Fig. 2.2 exhibits the inllucncc of an applied stress field on 
the hysteresis loo]) (Ilels. [10] and [11]). It is shown that a parallel compressive stress 
decreases the spontaneous and remanent polarization as well as the coercive field, but a 
conijjressive stress perpendicular to the field only slightly increases the coercive field. The 
temperature is also an important factor which influences the hysteresis loop. Figure 2.3 
shows that the remanent i)olarization Pr and coercive field Eo are functions of tempera- 
ture. It is shown that the whole hysteresis loop shrinks with increasing temperature. It is 
thus revealed that the hysteresis loop is a manifestation of the coupled electro-mechanical 
behavior of piezoceramics. 

In addition to the P-E loop, i)olarization switching leads to strain-electric field hys- 
teresis. A typical strain-field respon.se curve is shown in Fig. 2.4. The shape resembles 
that of a butterfly, and it is often referred to as the “btdierjiy loop". As the field is 
increased, the strain is no longer linear with the fiehl as domains start switching (Ref. 

I'.4|). 

The maximum positive fuild is limited by the breakdown of the dielectric, which in 
this case is the ceramic itself, and the v'alue of this Iwiiliruj field is 10 kV/cm, whereas 
the maximum negative field is limited by the piezoceramic depoling. When exposed to 
a high electric, field op])osito to the poling direction, the i>iczoceramic loses most of its 
piezoelectric capability, accompanied l)y large dielectric losses and poor efficiency. This 
is known as depoling of the f)iczoceramic and is accotnpanied by a permanent change in 
dimensions of the specimen. This is probably due to largo scale domain switching in the 
material. The DC depoling field of PZT-5H is approximotely 5.5 kV/cm (Ref. [8]). 

2.1 Constitutive Relations and Material Property 

/\.ri understanding of piezoelectricity begins with the structure of the material. To better 
explain it, let us consider a single crystallite from a }>olycrystalline ceramic. This crystal 
is made up of negatively and i)ositively charged atoms occupying specific positions in a. 
repeating unit or cell. The specific symmetry of the unit cell determines the possibility 
of piezoelectricity in the crystal. All crystals can be divided into 32 classes or point 
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groups (from 7 basic crystal systems: tricliiiic, monoclijiic, orthorliombic, tetragonal, 
rhombohedral, hexagonal ajid cubic). Of the 32 classes, 21 do not possess a center of 
symmetry, and 20 are j)iezoeiectric (one class, although lacking a center of symmetry, is 
not piezoelectric because of the combination of other symmetry elements). The lack of 
center of symmetry means that a net movement of the positive and negative ions with 
respect to each other as a result of stress produces an electric dipole. The ceramic, being 
composed of the random orientation of tlicse piezoelectric crystallites, is inactive, i.e., the 
effects from the individual crystals cancel each other and no discernible piezoelectricity 
is jnesent. A volume of the niaterial in which the ])olar moments of the unit cell are 
all oriented in the same direction is called a domain. The unpoled ceramic consists of 
many randomly oriented domains and thus Ims no net polarization. Application of a 
high electric field has the effect of aligning most of the domains in such a way that the 
polar axes of (.he unit cells are oriented as clo.sely parallel to the applied field as . allowed 
by the crystal strucl.ure. d’his i)rocess is called poHiuj a.nd it imparts a permanent net 
polarization to the ceramic. Since the polar axis of the unit cell is longer than the ol;her 
two axes, the reorientation of the domains also creates a permanent mechanical distortion 
Fig. 2.5 . The polar axes can be reoriented in cerlain specific directions governed by 
(,he crystal symmetry in the presence of a high electric field or mechanical stress. This 
process is called swUcInng. In the case of tetragonal PZ'Fs , either 90 degree or 180 
degree switching can occur (Fig. 2.6). The 90 degree switching is always accompanied 
by mechanical disl.orl.ion due lo the longer length of (.he unit cell along the j)olar axes, 
while the ISO degree swii,chlng can occur without any distortion. It is also known that an 
electric field can cause both 90 degree and 180 degree domain switching but a mechanical 
stress can only cause 90 degree switching. Conseciucntly, a mechanical stress can dcpolc 
a piezoelectric ceramic but cannot rei)ole it. The poled ceramics exhibits the direct and 
converse piezoelectric effects. 

2 . 1.1 Constitutive Relations 

For mechanical problems a constitutive equation describes how a material strains when 
it is stressed, or vice-versa. Constitutive equations exist also for electrical problems; they 
describe how charge moves in a (dielectric) material when i(. is subjected to a voltage, or 
vice-versa. Engineers are generally familiar with the mosl. common mechanical constitu- 
tive' eeiuations that ai)plir« for (weryda.y metals and pla.s(,ics. 'I’liis (Miuation is known as 
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Hooke’s law and is written as 


e = s • cr 

In words, this o(iuation states: Strain = Compliance * Stress. 

However, since piezoelectric materials arc concerned with electrical properties too, we 
must also consider the constitutive equation for common dielectrics: 

D = e-E 

In words, this equation states: Charge density — Permittivity * Electric field. 

Coupled equation: 

Piezoelectric materials comirine these two seemingly dissimilar constitutive equations 
iiitf) one coupled c((ua(,ion, written as 

c = Se ■ cr + d' ■ E 

E — d ' o d" £(r • E 

The matrix d contains tlic |)iezoclcc.tric cocflicionts for the material, and it appears 
twice in the constitutive e(|uation (the supetscript t stands for matrix transpose). 

In order to describe or model i)iezoclectric materia, Is, one must have knowledge about 
the materials mechanical properties (compliance or stilfness), its electrical properties 
(l)crmittivity), and its piezoelectric coupling proimrtics. The four state variables (e, a, D, 
and E) can be reai ranged l,o give an. additional three forms for a piezoelectric constitutive 
e(|iiation. Instead of the coupling matrix d, they contain the coupling matrices e, g, or 
q. It is possible to transform piezo constitutive data in one form to another form.- The 
four possible forms for piezoelectric constitutive eciuations arc shown below. The names 
of each of the form is arbitrary; they wore taken from the two dependent variables on the 
left-hand-side of each ec[uation. Note that the voltage and electric field are related as. 



Strain-Charge form; 
e. = Se ■ cr + d! ■ E 
JD — d ■ a -\- Sn ■ E 

Strain- Voltage form: 
e = So ■ (7 -h f/ ■ D 
E := -g ■ a + 6-^ ■ D 
Matrix transformations for converting 
into another are shown below. 

Strain charge to Stress charge 
Ce = 

e = d-sl^ 
e, = e^-d- si' • 

Stress charge to Sl.ress Voltage 
CD = Ce + e' • el' ■ e 
q = ■ e 


Stress-Charge form: 
o = Ce ■ e — E' ■ E 
D = e.-f. + e,-E 

Stress- Voltage form: 

a = Co ■ c - q'- ■ D 

E = —q ■ e A- ■ E 

piezoelectric constitutive data from one form 

Strain charge to Strain Voltage 
So = se ~ dEel^ -d 
g = el' -d 

train Voltage to Stress Voltage 
CD = sj/ 

1 = (J- Sp 

el' = el' -I- g ■ sj/ ■ .9* 


The subscripts in piezoelectric constitutive ecpiations have very important meanings. 
They describe the conditions under which the material i)roperty data was measured. For 
e.g., the subscript F on the comi)Uance matrix se means that the compliance data was 
measured under at least a constant, and preferably a zero, electric field. Likewise, the 
subscript a on the permittivity matrix e^ means that the permittivity data was measured 
under at least a constant, and preferably a zero, sl.ress field. 

A description of all mai.rix variables used in the i)iczoelectric constitutive equations 
is given as follows: 
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Symbol 

Object type 

Size 

Units 

Meaning 

a 

vector 

6*1 

N/m'-^ 

Stress comp)onents (e.g. <7i) 

e 

vector 

6*1 

m/m 

Strain components (e.g. ej) 

E 

vector 

3*1 

N/C 

Electric field components 

D 

vector 

3*1 

c/uT^ 

Electric charge density displacement components 

s 

matrix 

6*6 

m‘'^/N 

Compliance coefficients 

c 

matrix 

6*6 

N/m^ 

Stiffness coefficients 

s 

matrix 

3*3 

F/m 

. Electric permittivity 

d 

matrix 

3*6 

C/N 

Piezoelectric coupling coefficient for strain charge fori 

e 

matrix 

3*6 

C/m2 

Piezoelectric coupling coefficient for stress charge fori 

g 

matrix 

3*6 

n/^/N 

Piezoelectric coupling coefficient for strain voltage for 

q 

matrix 

3*6 

N/C 

Piezoelectric cou|)ling coefficient for stress voltage for 


The modified electio-thermo-elaslic constitutive equation including tlie nonlinearity 
due to polarization-electric field (P-E) interaction is given as (Ref. [7]). 

^ij ^ijk^k ^ f' ijO ( 2 . 2 ) 

Di = eijkfjk + SijEj -f (liO (2.3) 

7j = T (2.4) 

where, 

P— ^-Polarization 

The linear part of the stress is given as alj = Cijkitki — GijkEk — kijO, and the nonlinear 
part is = —PiEj 

2.2 Property Matrix for a poled Piezoelectric Ceramic 

A jriezoelectric cei amic has only oiie tyi)e of piezoelectric matrix regardless of the sym- 
metry of the constituent crystals. The ceramic is initially isotropic. This isotropy is 
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destroyed in tlie poling direction. In the direction perpendicular to the poling direction, 
the material is transversely isotropic. The elastic, dielectric and piezoelectric matrices 
for the cylindrical symmetry of poled PZT are shown below 



C|2 


0 

0 

0 


0 

0 

-ear 

C 12 


<^13 

0 

0 

0 


0 

0 

~^3l 

Cl3 

0 

Q3 

0 

0 

0 

e = 

0 

0 

— 633 

0 

0 

C/(/| 

0 

0 


0 

-ei5 

0 

0 

0 

0 

0 

C 44 

0. 


~ei5 

0 

0 

0 

0 

0 

0 

0 

CG6 


0 

0 

0 


e — 


Ex 0 0 

0 £i 0 

0 0 e.-, 


The subscrip(.s .1, 2, 3, 4, 5, and 6 refer to x, y, z, yz, xz, and xy respectively. 


2.3 Assumptions 

Before deriving tlic eriuations necessary to generate the hysteresis and butterfly loops, 
let us consider the assumi)ti()ns made in the. formulation. The assumptions made arc as 
follows (Fig. 2.7); 

1. The planar surfaces of the disk are stress free and in view of the thinness of this 
disk we may assume that the stress o^jindcy are zero throughout the specimen for 
all times. 

2. We further restrict our considerations to isothermal conditions so that tlie absolute 
temperature is constant and need not be included. 

3. Shear strains (75,2 — 7,:j/ = ^zx = 0) and shear stresses are assumed to be zero. 

4. Normal stress in z direction is Jtssumed to bo zero i.e., — 0. 

5. Piezoelectric material is transversely isotropic. 
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6. Electric field acts only in z-directioii i.c. Ex=Ey = 0. 

7. The [)oliiig direction of the piezo patches is along z-direction (see Fig. 2.7). 


2.4 Modt'is for hysteresis and butterfly loops 

In the present study two models have been analysed. These are: 

1. Complex Stiffness Model. 

2. Improved Domain Switching Model. 


2.4.1 Complex StifTiiess Model 

Imposing the necessary conditions (inentioned in the assmnption.s) in the constitutive 
Eejs. 2. 2-2. 3 , we obtain a set, of reduced constitutive equations given by 




-c*E 


-I- c*e. 


PzE. 


(2,5) 


D, — £*E, -I- e*f , 


Where, 

(I — 


e* = ... 

Cn4T!i2 


£ ' — £3 


Op2 

Fn iT’i-;, 


But we know from Eq. 2.1 (Iia(,, 


D, = P, -I- £0 • E, 


= 4 » / 2 = Dz — £0 ■ Pz 

Substituting 0% from E([. 2.G into Eq. 2.8, 


( 2 . 0 ) 


(2.7) 


( 2 . 8 ) 
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Pz — (e* — Eo)Pz -|- C*c^ (2-9) 

Thereforo, sul>st,il,uuig r'oni Eq. 2.9 inl'.o Eq. 2.-5 and putting < 7 ^ = 0 we get the 
following relation for strain 

and Eq. 2.9 can be written as, 

Pz ~ e\Ez -f- E*e.z (2.11) 

where, e* = <r‘ -- Eq. 

If we assume Ez — EoP"’'', wlicre ui corresponds to the excitation frequency, the strain 
variation exhibits a, 2uj variation due to tlie E^ term. 

Strain in tlie above ef|nation is a double harmonic function of electric field and this 
nonlinearity is resironsible for the generation of butterfly loop. The material properties 
of PZT-5fI are given in Table No. 2.1. Using these values, the Eqs. 2.10 and 2.11 can be 
written as 

_ 151G8 * 10"''^/';^ * ^.{2w/. -Q.(mq gt(^r-o.o.32,';) 

N- -- T:,M9 + 10'" - i9.9/io 

P, 1 5 1 08 * 1 0 Eo * -i- i9.9To * (2,12) 

The factor ~ 10108* 10 is found to be negligible compared to e* ~ 19.9 (multiplying 
factor of the linear term) and hence the butterfly !oo{) can only be generated at higher 
electric fields. I(. was found in the literature survey that butterfly loop occurs at even 
low electric fields. So, we (h'.vclop a model which would generate both the hysteresis and 
butterfly loops at even low electric fields. 

2.4.2 Improved Domain Switching Model 

/ 

This model incorporates the conccirt of domain srniic.hmg in the complex stiffness model. 
It is a simplified ve.rsiou of the model developed in Ref. |14|. It is generally agreed that 
the hysteresis loop is a cons(xiuei\ce of the dela.y in resironse of domain switching. From ' 
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a microscopic viewpoint, the niechaiiisni of doinaiii switching are not well understood. 
Hence it is dillicnlt to formulate explicit theories based on the present knowledge of 
tliese mechanisms. On the other hand, polynomial relations which partially match the 
graphical rcj)resen(.ation ol the loop have not proven to be completely satisfactory for one 
reason or the another. For instance, they yield mniti-valued relations and inadmissible 
states between the electric displacement and the electric field, and they do not explicitly 
account for the time dependent effects of domain switching. 

It is our intent to offer a macroscopic theoretical basis for the existence of the hysteresis 
loo{). This theory is motivai.ed by the fact that the effective number of dipoles aligned 
in the direction of the applied field is altered as domains switch, and this number obeys 
a rate law which dei)ends linearly on the electric field. It was observed that the specific 
nature of the loop depends not only on the amplitude and the i)eriod of the cyclic field, but 
also on its entire temi)oral history. These considerations point to the fact that time must 
be an important factor in any theoretical consideration of the loop and this explains the 
inclusion of the rate law (in combination with the constitutive equations) in the modeling. 

Therefore the modified electro-clastic constitutive equations after having included the 
imiiortant parameters i.e. the rate law and number of di])oles is given by 

(rij = cijkifki - CijkNEk - PiNEj (2.13) 

Di=^(kjkNejkPeijEj ' (2.14) 

+ (2.15) 

Where, N is the effective number of aligned dipoles and it obeys a rate law given by 

N-\-aN^bE (2.16) 

Where, a and 6 are constants. 

Let, E = E^^E^ 

and N = iVe“'- 

substituting (.he above twf) lorms in rate law i.e. E(i. 2.1G we get 
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(2.17) 


J\[ =z .. g’-H—/’) 

v 


Where, b = hE{) 


(j) = t!\n ’(— ) (2.18) 

a 

Substituting and b=0.1 in Rq. 2.17 the expression for N can be written as 

TV = 

sj 1 -I- oj'^ 

Let A — — /f/ .j then, 

N = (2.19) 

The reduced set of constitutive equations after imi)osing the necessary conditions is 
giv^en as 

.r, =. -e*E,N - I\ E, N (2.20) 

D, = 6*7?, -t-cT,At (2.21) 

From Efjs. 2.7 and 2.21 p(,>lari7,a,tion can be written as, 

/ 2 = (e* — i^E}Ez + E CzN 


Pz = e\Ez + EczN ' ( 2 . 22 ) 

Substituting IV, bom Eq. 2.22 and imposing - 0 {stress free condition) we obtain 
strain in terms of N and electric field. The relation is given as 


e^E^zN + EEzN 

c* — eAEzN"^ 


(2.23) 


Substitutijig the values of e\, E, E, and expression for N (Ec]. 2.19) the Eqs. 2.23 
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and 2.22 can i)C written as 




AP^l t 15108 10 * ,,7{.3to(-n.o,3.i(>-^) ^ ^ ^i(2^t-om25-4>) 

5"3 1 irriO"' * 19.9 * 


(2.24) 


= 151G8 * 10-'^ * jg g ^ ^ ^ g«M-o.o325-^) 


Strain is nonlinear in terms of electric field (Eq. 2.24) and this nonlinearity is respon- 
sible for the i>rescnce of butterfly loop even at low electric fields. 


2.5 Results and Discussions 

Hysteresis (j)olariza,tion-el(M:ti ic fi('ld) and butterfly (strain-electric field) loops are plotted 
for the two different models discn.s.sed. The variation in nature or shape of these curves 
with change in different parameters is studied through the„sc plots. Figures 2.8 and 2.9 
represeni. the plol.s for complex stitfiKi-ss model. Plots 2.8 (a), (b), and (c) show the 
variati(,)n of P-1'1 and r — ft' loops with change, in freepiency. Electric field is fixed to a 
value of Eq-- Io HYui in tlu'se three plots. It is observed that there is no variation in 
these loops with variation in fixupiency. In Figs. 2.9 (a), (b), and (c), frequency of the 
system is kept, constant (i.e., J —■ 0.01 llv) and the effect of electric field on the loops 
is studied. Diastic changes are olrserved in the loops with variation in the magnitude 
of electric field. The strain-electric field plot transforms into a double loop termed as 
“huUerJly loop" ni higher electric fields (of order of IO''^V/m). Even the hysteresis (P-E) 
loop transforms in a. similar manner. 

Figures 2.10-2.13 reincsent I, he behavior of hysteresis aiul butterfly loops with change 
in parameters ch'fining the domain switching model. They clearly depict the effect of 
frcciuency, magnitude of electric field, parameters a and 1) defining the rate law, on the 
P-B and e - E looi)s. Froni Figs. 2.10 and 2.11 it is observed that the shapes of butterfly 
and hysteresis loops change drastically with change in frequency but there is no change 
observed with variation in electric field (from KPV/m to 10® V/m). From Figs. 2.10 
(c) and 2.1 1 (c) it is (.)bservcd that higher frequencies the. butterfly loop transforms to a 
loop similar to tha,(, obtaiiuxi for electrostrictive materials (ll.cf. [3]). In developing those 
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developing these figures (Figs. 2.10-2.11), the values of a--l and b=0.1 are used. The 
effect of parameter “6” on tlie loops can be clearly studied from Figs. 2.12 (a), (b), and 
(c). The butterfly loop cx])ands with increase in b and the hysteresis loop transforms 
into a triple loop at higher value of b (for b=0.5). From Figs. 2.13 (a), (b), and (c) it 
is ob, served that the buttcrlly loop shrinks gradually with increase in “a” and not much 
change is observed in the hysteresis loop.' 


2.6 Summary 

A niathcmatical model has been dcvelojrcd to represent the interaction of P-E (polariza- 
tion and electric field) and domain switching has been developed to obtain the hysteresis 
and butterfly loops for a piezoelectric (PZT-5II) material. Basically two models are 
studied, one with domain swi(, clung and the other without it. The effect of inclusion 
of domain switching is shown to represent the physical phenomenon observed for PZT- 
niatcrials. The effect of varying the parameters in the rate law (frequency, electric field, 
a and b) on the P- E and c — E loops is analysed. 
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Figure 2.2: Change in hysteresis curve with respect to parallel compreSvSive stress 
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Figure 2.5: Dimensional changes of a poled wafer with an api)lied voltage. The dimen- 
sional changes are greatlj'’ exaggerated for clarity. 
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Figure 2.G; Switching 
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(a) f = 0.01 Hz 



Electric field (V/m) 

(b) f == 10 Hz 








Strain (m/m) ^ ^ Strain (m/m) 


Electric field (V/m) 


• — • — I 

JO® Electric field (V/m) ^ iq® 

(a) E = 10^ V/m 






(c) E ICP V/m 


Figure 2.9: Variation of strain, polarization, and electric field with change in magnitude 
of electric field (f-d),01 Hz) (Complex stiffness model) 
























Figure 2.12: Variation in hysteresis and butterfly loops with change in parameter b 
(F-10*^^ V/m; fro(iuency- 0.01 Hz ; a— 1 ) (domain switching model) 









(c) a=2 

Figure 2.13: Variation iii hysteresis and butterfly loops with change in parameter a 
(E"10^ V/m; frequency=0.01 Hz ; 1)=0.1 ) (domain switching model) 
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Chapter 3 


HEAT TRANSFER ANALYSIS OF 
SMART ACTUATOR 


Heating of smart actuators is a major problem in tlx; industry. We, in this study have 
tried to analyze tlie heat transfer that happens within tlic actuator i.e., as to what is the 
heat generated within the actuator and how it gets dissipated and what is the maximum 
temperature reached when conduction and convection from cjkI faces; conduction and 
convection from all faces arc considered? These are some very important questions that 
have been addre.s.sed in this study. 

The hysteresis effect represents the energy lo.ss. This energy loss is responsible for 
lieat generation within the piezo material. The heat so generated, leads to heating up 
of the piezo layer tremendously, .sometimes causing .serious damages to the actuator. So, 
measures have to be taken to restrict this energy lo.ss to a mininmm. In this formulation, 
the temperature profdc within the actuator is found out. For this the ba.sic taws of heat 
l.ransfer are considered (Refs. [25|-[27|). 


3.1 Modes of Heat Transfer 

The basic heat transfer (or heat di.s.sipatio]i) phenomena considered here are namely: 

i. Conduction : When there exists a temperature gradient within a body, heat energy 
will flow from the region of high temperature to the region of low temperature. This 
j)henomenon is known as conduction heat transfer, and is de.scribed by Fourier’s Law 


53 



(named after the French i)h 3 \sicist Joseph Fourier), 

q = -fcvF 

'I'his equation detertiiines the heat flux vector q (fjuantity of heat/unit area) for a 
given temperature profile T and thermal conductivity k. The minus sign ensures 
that heat flows down the temperature gradient. For problems where the tempera- 
ture variation is only 1-dimensionaI (say, along the /.-coordinate direction), Fourier’s 
Law of heat conduction simi)lifies to the scalar equations, 

q = -k^ Q^-kA~ (3.1) 

dz oz 

In Uie above equation on the I’igbt, Q represents the heat flow through a defined 
(uoss-seclJonal area A, and it is measured in Watts, 

Q — f q ' (lA 

Ja 

The ID heat flow eciuation through a material thickness L is given as, 

Q = M(7’_'r2) (3.2) 

Jj 

where T| and T 2 arc the temperatures at the two boundaries (i.e., at z=0 and 
7,-f.). 

2. Cantiection : Heat energy transfers between a solid and a fluid when there is a tem- 
perature difrerence between the fluid and the solid. This is known as "convection 
heat transfer". Generally, convection heat transfer can not be ignored when there 
is a significant fluid motion around the solid. The temperature of the solid, due to 
an external field such as fluid buoyancy, can induce a fluid motion. This is known 
as "natural convection" and i(, is a strong function of the temperature difference 
lietween the solid and the fluid. Blowing air over the solid by using external devices 
such as fans and pumps can also generate a fluid motion. This is known as "forced 
convection". Heat transfer due to free convection is described by Newton’s Law of 
Cooling, 


Q = /,,/}(7;„ - T’„) = hA ■ AT 


(3.3) 
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Thf! rale of hnal transfenod Q to the siurouiidiiig fluid is proportional to the object’s 
exposed area A, and the difference between the object temperature and the 
fluid free-streani temperature 7jj. The constant of proportionality “h”\s termed the 
convection hcat-traiisfer coeflicicnt. 


3.2 Analysis of Heat lYansfer within the Actuator 


The smart actuator is considered in actuation mode i.e., oidy external voltage is applied 
on the actuator (Fig. 3.1). The actuator consists of alternate layers of metal-piezo 
arrangement. The amount of heat generated is calculated from the area of hysteresis 
loop. 

The heat generated/unit cycle/unit volume, ij is given as 


</ = 




(It 


( 3 . 4 ) 


where P is the polarization which is given by E(|. 2.12, Tp is the time period, and .L 
(--/fo sin wi) is the electric field. Here Eq is taken to be 10^’ V/m. 

Two distinct formulations are described in the following section. 


3.2.1 Mode of Heat dissipation : Conduction and Convection 
from end faces 


A layer-by-layer fortnulation is presented over here. 

dlie general 3-D steady state heat couducl.ion e(iuatic>n (obtained from the conseiva- 
tion of energy principle) is 


d'’-T d^r ^<1 

Therefore, Hie 1-D steady state heat conduction equation is given as 


c?'^r q 

dz^ ^ k 


= 0 


r 

where (7= energy generated per unit volume, W/nC 


k “ Thermal conductivity 


( 3 . 5 ) 


( 3 . 6 ) 



The smart ac(.uator is a stack of alternate layers of Jiictal and piezo material. Here, 
we consider the SJiiart actuator to l.)e a stack of "n” layers. 

1. First Metallic layer (Fig. 3.2) 

Since there is no heat generation in tliis la.yer, the governing differential equation is 


d'^T 

dz^ 


= 0 


(3.7) 


subjected to boundary conditions that on face z=La/ the metalic layer is subjected 
to temperature T 2 and on face z=0 it is open to the atmosphere for convection (Fig. 
3.2), 

i.e., at z=(), = h{T - Ta) and at T = T 2 

where /ca/ and arc the thermal conductivity and length of the metallic layer 
Solving E(i. 3.7 and using the boundary conditions, we get the following rela.tion 
for tcni])ora.ture profile 

T{z) = CyZ d- C2 


where, cq 




2. Last metallic layer (n th layer) 

For this element, the l)oundary conditions are: 

at 7-Lm, -k^f- = h{T - Ta) and at z=0, T = '1\ (Fig. 3.3) 

where Hm and are iho iJiennal conductivily and lengtli of the metallic layer. 

The tenii)erature profile in this cavse is given as 


T{z)=.c,z~\-C2 


where, c\ 


JjtTzJj 

//i) 


and €'2 = 


3 . Metallic layer (between two j.'iezo layers) 

The boundary conditions in this case are 
at z-0, T = 7’, and at z^Lm, T = Tz (Fig. 3.4) 
Then the tcmi)eraturc profile is 


T{z) = C\Z c'z 


where, rq — 




ami C2 = 7 j 


5C 



4: Piozo layer 

The differential equation in this case is slightly modified (as compared to the other 
two) to include the internal luiat generation term and is given as 


(Tr 

Oz^ 



The boundary conditions in this case are at z=0, T = 2\ and at z=Lp, T = T2 
(Fig. 3.5) 

where kp and Lp arc the thermal conductivity and length of the piezo layer. The 
temperature profile is obtained as 


T{z) = r, + {T2 - 1\)— + 

i-J-n 


qq, 

2kp 


- - 

p 


When a comlunation of MctaI-Pi(!zo-Mcta,l layers is taJeen for (;xami)le, (MPM block as 
shown in Fig. 3.G) then the temperature profile of the entire bar is given as a function 
of temperatures and (i.o., te?ni)cratures at second and third junctions as shown in 
Fig. 3.6). The heal, continuity conditions at the iuterfa.ccs are respectively 


-k 


Af 


dz 



-kp 


dT 

7z 


{at z = 0)i 


{Metal — Piezointerface) 


-kp— (a/: z = L„)i — —k^ -f- {at z = 0),-(.i {Piezo — Metal interface) 
dz dz 

The subscrii)t denotes the layer number. 

Similarly any number of alternate layers of metal and i)iezo can be taken for analysis. 
In i,he present formulation a bar of 25 layers (made up of alternate layers of metal and 
piezo i.c., MPM. ...PM) is analyzed. 

3.2.2 Mode of Heat Dissipation : Conduction and Convection 
from all faces 

In the earlier formulation convection from the cylindrical surface has been neglected. The 
|)rocedure for inclusion of convection from all the faces is described in this section. 



VVf! approach I, ho problem l)y making an energy balance on an element of thickness 
dz. Thus 

I'aiergy flow t.luongli I('ft la.c(i -- energy flow (.hrongli light face d- energy lost by 
conveci ion (eiungy balance .for a metallic hvyer I'ig. 3.7) 

Qz — Qz+dz + }lPMdz{T — Ta) 

where, Pm is the perimeter, = -k^A^ and Qz+dz = -kMA{^ + ^dz) 

And energy balance for piezo layer (Fig. 3.8), 

Qz -f- Q, = Qz+dz + hPpdz{r - Ta) 


whore Qg(h.cal. (/r+icraicd) = qAdz, Pp is the perimeter, Qz = —kpA^ and Qz+dz = 

-1^ S'-fc) 

Substituting for Qz and Qz^dz erie obtains 

dz‘ kA ' ■ 

The layer- by- layer fornndation is given as 
I . First metallic layer 

The govei ning differfMitial equation applicable l.o this case is given exs 

on' iiPm „ 

Dz'^- IzmA^ 

Let 0 = {T - :/’„) and m.l, = 

'bhen it follows that 


d^O 2 n 
0 ^ - - 


(3.8) 


The boundary conditions arc given as 

at z;--0, — hO and 2 = Lm, 0 = O2 (Fig. 3.2) 

Solving cq. 3.8 and using the above given boundary conditions, wc obtain a tem- 
perature ixrofile as follows 


0{z) = 


0< 


2 


co.sh(7n/ifLAf j I- 


(X)sh(?/?.fl/2:) -b 


h 


iriMk-M 


■ sinh(777.A./i) 



2. n th metallic layer (Fig. 3.3) 

For this element the bouiKiary conditions are given as 

The temperatnre profile for this case is given as 

^ - z)) + - z))] 

‘ cosh(rnA^LA./) T sinh(7nAfLM) 


3. Metallic Ia,ycr (between two piezo la.ycrs) 

The boundary conditions in this case are 
at z=0, <? = (?[ and at z—Lm, 0 = 02 (Fig. 3.4) 

Then the temperature profde is 

+ 01 ■sinh( 7 n,Af(TAr - z)) 


4. Piezo layer 

The (lifTcrential eqnai.ion in this case is slightly modified (as compared to the other 
two) to include the internal heat generation term and is given as 


Oz'^ 


77?,pd + 


up 


= 0 


where, 77!,|, 1-^. 

'Fhe boundary conditions in this case are at z=0, 0-0] and at z-Lp, 0 = 02 (Fig. 
3. ,3). Hence tin; tom])eiaturc profile is given as 

{01- i-^-) smh{mriLp - z)) + {O2 - - ^,t) smhjrnpz) ^ g 
^ sinh(m.pL7.j kpm], 

Whe.i a combination of these bayers is taken (for example MPM block as shown in Fig. 
3.6), then the temperature {)rofile can be found out in terms of (?2 {T 2 —Ta) and 0z (T3— Ta). 
These quantities are obtained by imi)osing the following heat continuity conditions 

-A;a/~ (ai 2 : = LM)i-i — ~kp— (ai ^ = 0)* 
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{niz =-- 0)m.| 

Similarly any mimbcr of ;vl (.ornate layers of metal and piezo materials can be analyzed. 
In the present foi nmlation a bar of 25 layers (made np of alternate layers of mel.al and 
piezo i.e., is a.naly7,ed. 

3.3 Results and Discussion 

Using the above heat transfer formulation, the temperature profile of a piezo stack actu- 
ator acted on by a harmonic, elecl.ric field (/?o~fO^V/in) is studied. In the present formu- 
lation a bar of 25 layers (made u{) of altcrna(,e layers of metal and piezo i.e., MPM....PM) 
is analyzed. I’ln; metallic' layers are of length 0.5 mm and the piezo layers arc 1.5 mm 
thick. The material and geometric data are given in 'fables 3.1 and 3.2. 

The tem})eratine profile within the actuator for various cases is shown in Figs. 3.9 
and 3.10. The results can be summarized as follows 

1. 'Fhe temperature luofile obtained in all cases is .symmetric about the center of the 
l,ar (Fig. 3.9). 

2. for (,lie ease, where convection from the surface area of t.he bar is neglected, it is 
observed (.lia(, for an inc.rea.se in fre.<iueuc.y by 10 llz the maximum temperature 
reached by tlu', l>ar goes up by ~28"K as shown in Fig. 3.9. 

3. For the case, where convecl.ion is allowed to happen from the surface area of the 
bar, t hen it. is observed that for an iucrcfisc in fref|uency by 10 Hz the maximum 
temperature reached i)y the bar goes uj) by ~2'’K as shown in Fig. 3.9. These results 
indicate (.hat inclusion of convection (from the surface area of the bar) reduces the 
temi)oraiure increase sul)stantially. 

The same calculal.ion.s vvere repeated for the case where .steel is chosen as metallic layer. 
The temperature profiles obtained were found to bo. similar as in the earlier case (Fig. 

3.10). 
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3.4 Summary 


Heating of pio7,o ac.tnatoi.s is a major problem in the industry. An attempt is made 
to fornndate f.lio heat transfer |)roblcni applicable for piezo actuators. Through this 
formulation we come to know the maximum tcm])erature reached by the actuator. The 
effect of inclusion of convection mode of heat dissipation in addition to conduction is also 
brought out. It was observed that for the same frequency the maximum temperature 
reached was much higher in the former model (model excluding the convection from the 
surface area of the bar) as compared to that obtained by latter model (model including 
the convection from all faces of the bar). F'rom the results it was found that the profile 
was symmetric about the cent, re of the bar. 
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lal)lo 3.1: Material Properties for PZT-5H 


GPa 

'T1 

<U2 


C33 

C 44 

^'00 

138.33 

93.10 

95.00 

131.26 

21.05 

21.05 


C/m^ 

^31 

T33 

Cl 5 

-4.0329- 1-0. 13131 

lG.5815-0.539i 

12. 29-0. 4i 


F/ m 


£3 

(15293-510i)*Hr^''^ 

(15028-501i)*10'^^ 


Table 3.2: Data u.secl for heat trairsfer analysis 

Diameter of the bar = f mm 
Length of the bar - 24.5 mm 
Number of layers = 25 (alternate MP..M) 

Thickness of the metal layer — 0.5 mm 
Thickness of the picz.o la. 3 'cr := 1.5 mm 
Ambient tempcratiire = 300 K 
Thermal conductivity of A1 237 W/mK 
''I'hermal conductivity of Steel = 23 VV/mK 
Thermal conductivity of PZT-5II = 1.8 W/mK 
Heal, transfer cocflicicnt for natural convection by air = 5 W/m K 
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I'igurc 3.4: Metallic Ia.yer (between 2 piezo layers) 
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Figure 3.5: Piezo la.ycr 
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Figure 3.G: MPM block 
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Chapter 4 


MODELING OF SMART ACTUATOR 


In Uiis chapter a finite element formulation is carried out to obtain a relation between 
deformation, potential and extcrtial load acting at the ends of the smart actuator. The 
smart actuator is a circular bar which consists of a stack of alternate layers of metal 
and j)iezoeIectri<; material. The smart bar acts Jike a sensor (Fig. 4.1a) when induced 
potent, ial is measured due to deformation caused by the external load. Similarly the bar 
can be used a»s an actuator (Fig. 4.1b) when external potential is applied across piezo 
patches to create deformation and force at the cjids of the smart bar. 

PZT-5II is (,hc material chosen for the piezo layer in the present formulation. Bar 
model is shown in Fig. 4.2. The assumptions considered in this formulation are : 

1 . The planar surfaces of ( he disk are stress free and in view of the thinness of this 
disk we ma,}- assume that the s(,rc.ss ay-jinday arc zero throughout the specimen for 
all times. 

2. We further restrict our considerations to i.sothermal conditions .so that the absolute 
(.emixiral.urc; is constant a, ml need not be includetl. 

.3. Shear st.ralns (7^^ — 7.,;,/ = 'Jz-x ~ 0) a»d slicar stresses arc assumed to be zero. 

4. Piezoelectric material is transversely sotropic. 

5. Electric field ad,s only in z-diredaon i.e. Er= Jf, = 0. 
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Finite Element Formnlation:- 


Tlic finite eleiiH'iil. fonniilatioii is based on using the variational principle. The variational 
form of the liiu'at electro-c'lastic prolilcni.is given as follows (Ref. [6|) 

ajjSfijdV - ^ DiSEidV = l^fMidS 

For a piezo layer, in the absence of traction force, the variational formulation is given 
as 

o’;Se,dV - VJE.dV = 0 (4.1) 

In this stud}^ wo ha.v(^ not considcroc! the noidiiicarity due to P-B interaction. There- 
fore tfio roduood ('onstitutive ocmations for a linear ea.sc3 are given by, 


-c*E, F c*p. 

(4.2) 

: -1- 

(4,3) 


Where, 


r" rrr ^33 


2r 


13 ..,. 


f’ 1 t \ (' 1 7 . 




rj j I C12 


:i3 


e* = d ■ 

t')i 

Substituting the constitutive IRi 4.2 and 4.-3 in the variational formulation (Eq. 4.1) 

implies 


- e*E,)Sc,dV - /^(e*^z - e*c.,)5E,dV = 0 
Substituting for r, = and E, = Eq. 4.4 can lie written as. 


(4.4) 




£'~+e*^')^dl/ = 0 


f , (if)?/. . f f 

/,■*'■ Ih. jdTz id ^ Oz ' " Oz' dz 

Integrating over the cross-sectional area, we have 
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+ c’ 


d^^dSu 


)~;—dz -f- A 


Oz ' d. 



^du.d6^ 

&>-§r 


dz — 0 


A fiiute ol(Miu!ii(, fonmilal.ioit has been clcvcloi)f!d using a la,yor-by-layer model satisfying 
(lisi)laccmcnt continuity. Both axial displacement u and electric potential # are assumed 
to be quadratic functions of z. ddie element used for modeling the piezo patch has three 
nodes. Each node has two degrees of freedom. These correspond to one extensional 
displacement in z-direction and one electric potential in poling direction. This results in 
a stiffness matrix of the order of G*6 for the i)iezo patch element. Whereas the element 
for the non-piezo material ha,s three nodes, having one degree of freedom per node. These 
corresfjoud to one (vxtensiona.1 degree of freedom in the poling direction. This results in 
a stiffness matrix of the order of 3*3 for the non piezo material (Fig. 4.3). 

The axial displacement and electric potential are given in the discretized form as 

”(70 " EjL 1 »j dui = dUiNi 

dd>i^5<hNi 

where Uj and <T>y are nodal degrees of freedom (Lagrangian) and Nj are the shape 
functions. 

yV, ... (I .... „ g.) 

N. - (ft)” - 
ft - ( ■- ft) 

Substituting u a.nd d> in tlie virl.ual work formulation we get 


fi-r ONi * , 


= 0 


(4.5) 


On substitni.ion of u, <I>, (hi, fFI>, wo get the local ,stijjneas matrix for the piezo layer, 
which is given as, 
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Isl -8c* ^ 7^ -8e‘ e* 

8 3 3' 3 T 

nic-l ICc* -8c* -8e* ICe* -8c* 

3 3 3 3 3 ' ~3~ 

-8c* 7c* c^ -8c‘ 7e‘ 

3 3 3 3 3 3 

-8e* cl -7e* Be* -e* 

3 3 3 3 3 3 

-8c‘ ICc* -8e* Sc* -ICe* 8c* 

3 3 3 3 3 3 

Sl -8e* Tel z:£l 8el -7c* 

3 3 3 3 3 3 


Tlie nodal (iisplacnmonf, v(x:(-or is 


Those are also caJlcd as Primary Variables 

Element stiffness matrix for a metallic layer (Ref. |24|) 

7 -8 1 ’ 

-8 IG -8 
1 -8 7 

where If is the '\7)ung’s Modulus of Elasticity. 


'The nodal displaceuHUit vector in this case is givcni as 


In deriving the clement cciuations, wo isolated a typical clement from the mesh and 
fornndal.ed the variational problem and developed its linite element model. To solve 
the total problem, we jnit the element back into their original position. In doing this 
before discretizaXion, w(', imirose continuity of the primary variables at the connecting 
nodes between elements. Continuity of the primary variables refer here to the single 
valued nature of the solution. I'hns the assembly of elements is carried out by imposing 
continuity of primary variables at connecting nodes i.e., xif, = (the last nodal value 
of the element f> is same as the first nodal value of the adjacent clement e+l). After 
imposing these conditions we get a. (li.scretized linear electro-elastic oriuilibrimn equation 
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for the entire smart jiiezo aetnator which is given as 

|/l'l{a) = {F) 

vvlinrp |A'| is l.lis glolwil slilDioss miitrix, («) is vector ot global degrees of freedom, 

and (F) is the global load vector. 

A finif.e element code is written to carry out the formulation and the required results 
are generated for sensing and actuation cases. 


4.1 ResuKiS and Discussion 

! 4.1.1 Sensing inode 

! The results general.fvl from finite element formulation show that in sensing case the tip 
i donection and potential induced can bo expressed by the linear relations as 

! 

ii'iip = e^Frs (4.6) 

and 

dhriffuro/ = l>[Frs (4.7) 

where (p and />| ai(' con.stants. I'ln' external mechanical load Fps is taken to be positive 
when tensih' and negal.ive when it is compressive. Since both uup and ^ are linear 
functions of tin- sensing Jotcc, the constants Oi and /q arc found out by applying a unit 
force at i.he ends of the ba,r i.e., /')>.s = IN. For a given set of data (Table 4.1 and 4.2) 
the values of fp and hi are found to be 

fp -- 4.(1S28() + 10 ^ and /p ~ .'51. .582 (when Aluminium is chosen as the metallic 
layer). 

4.1.2 Actuation mode 

In actuation mode, (.he acttial.or is acted upon by external mechanical load Fpa {Fpa is 
taken to Ire positive when tensile and negative when it is compres,sive) and an external 
potential as slnnvn in iMg. d.lb. The tip deflection is of the actuator can be expressed as 
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’’hip — C\PrA + C2^'^rA 


( 4 . 8 ) 

wlicic (’i and '"O c.oiisl.aid.s. Ihc cous(,aiil, ci is obl.aiucd by applying a unit force 
and 7-cro volt, age (/'V.-t 1 /V and VVa = OV"") and while obtaining the constant C 2 it is 

assumed that /'Va ” ON a,nd = W. For the given set of data (Table 4.1 and 4.2) the 
values of C[ and are 

Cl = 5.5 * 10 ^ and c^ — —4.488 * 10“'’ (when Aliiininiuni is chosen as the metallic 
layer) 

In this case, piezoelectric ceramic actuators exhibit a decrease in their displacement 
for a given cxcil.aiion voltage as they are loaded, d’his relationship can be seen in the 
above ecjuat.icm, for t.he (uase where As the load is increased, the displacement 

eventually reacluis zero, for a given applied potential lA^i. This load, where the net tip 
deflection of the jiic'zo sta,(4< actual, or is zero, is termed as the block force. 

The restdts obl,ained from the finite element analysis are plotted as shown in Figs. 
4.4-4.10. The resnlls can be snmniarised as follows: 

1. Plots 4.4 a,b, and c .show the variation in deflection and induced electric potential 
wit.h change in sensing force applied at the ends of the bar. The deflection varies 
linearly along tin' s|)an of the piezo stack sensor and potential induced varies linearly 
within the pii'zo layer. 

2. In tlu' case of acl,nat,ion Plots 4.5 a, I), and c shows the variation of deflection 
along l,he bar for different a})plied potentials. Deflection (elongation) varies linearly 
wit.hin the pii'zo layer and remains constant within the metallic layer. 

3. Idot 4.5 represents tin' linear relationship between the applied potential and required 
blo(.:k force. It is observed that a potential of -200 V requires a block force of -1.62 
N .so as 1,0 get zero deflection at the tip. 

4. Pl(.)t 4.7 gives the nodal deflection when a potential of -123 V and the corresponding 
block force of -1 N is airplied. It is clearly seen that the jiiezo nodes (denoted by 
f,he ujqier point, s) experience elongation and the metallic nodes (denoted by all the 
lower point, s), e.xperience compression. The. net deflection at the tip is observed to 
be zero. 


73 



of the deflection in the piezo and metallic layers are found to be different. Figure 
4.9 show.s the variation of deflection along the bar when a voltage of -500 V is 
applied acro.s.s the piezo patches. Deflection varies linearly within the piezo layer 
and remains constant in the metallic layer. Figure 4.10 shows the relation between 
the potential applied and the rexpiired block force. 


4.2 Summary 

A smart actuaf.or/sensor consisting of alternate layers of mef.al and piezo material is mod- 
eled. A finite element formulation is undertaken to obtain a relation between deflection, 
])otentials and external mechanical loading for both sensing and actuation cases. 
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Table 4.1: MaI.erial Propcii.ics for PZT-5II 


GPa 


'T2 

-_^13__ 

^■33 

'■'44 

'TO . 

08. IG 

05.00 

131.20 

21.05 

21.05 


C/ny^ 

'’31 


1L5 . 

-4.0320 

ir)..58l5 

12.20 


r F/m 

ei 

f.i 

15203"''Hr^'^ 

15028*10'^''^ 


T<U)1(' 4.2: Da.i.a for uiocleliug t,lto actuator 

Dhaii^cirof tlic bar = 1 lum 
IjCiigtli of tlic bar “ 24. b initi 
Nintibcr of bayers used — 25 (alternate MP..M) 
Thickness of the metal la.yer = 0-5 mm 
Thickness of the piezo layer = 1.5 mm 


Young’s modulus for metal layer 
For steel E - 208 GPa 
For Aluminium E —'70 GPa 
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H Metal 

(a) Pi(''/,{) Sl.ack vSensof 




Metal 


\ 1 Piezo 


(b) Piezo Stack Actuator 


I'^igurc 'I.I: I'iczo Stack Sensor a.n(l yVctnator 



[•'igiire, 4.3; Finite element model of smart actuator 


7G 




■I coof (1 i n a t c ( Tn ) 7 , coord i na tc ( m) 


(a) Fpg= -I N 




0 

-5c-07 
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~4(> 06 
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•U 

6 -300 
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• 4 5c-06 

0 


..I 

0{>05 0 01 0.015 0,02 0,025 

7. cooidiriatc (tn) 


-350 



(c) Fpg--. -.10 N 


Fif’UK’ 4 . 4 : Vnriril.ioii of niul inducod pototil.ifil for different exteiniil mechanical 

load (Sensing ea.se inet.allic hiyer-ahuniniinn) 
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ignre 4.5; Vnrial.ioji of dofiocUon along tho span for different, applied potentials (Actu- 
doii case) 
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Block Force (N) 


A1 

Figure 4.G: Plot of f)o(,eiitial Vs block force 



AI 


Figure 4.7: Plot, of deliecUou along Uie bar when block force is applied 





Figure 4.8: Seu.sing mode (meUillio Ia.yer-sl,eel) 
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Potential (Volt) Detleciion (m) 



0.005 0.0 1 0.015 0.02 0.025 

7 . coordinate (rn) 


$ -500 V 


Figure 4.9: Plot for actuation case (metallic layer-steel) 



Olock Force (N) 

Steel 

Figure 4.10: Plot of potential Vs block force 
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Chapter 5 

VIBRATION CONTROL IN 
HELICOPTERS USING ACTIVE 
VIBRATION CONTROL (CLOSED 
LOOP POSITION FEEDBACK) 


The various sources of vibration in helicopters are the rotors, engine, and transmission 
system. Of these, the primary source of excitation is the main rotor system, which oper- 
ates in a complex aerodynamic environment. The rotor blade loads are generated by an 
interaction of the physical (luantities associated with the aerodynamic, inertial, and struc- 
tural proi)erties of the blade. These vibratory loads are transmitted to different parts of 
the fuselage through a complicated load path and cause discomfort to the pilot and crew, 
equipment deterioration, fatigue damage to the structure, and increased maintenance. 
The adverse effects of these vibratory loads increases with increase in forward speed of 
the helicopter ajid also the cumulative fatigue damage to the structure increases with 
higher utilization of the vehicle. As a result, these vibratory loads restrict the operation 
and efficiency of the vehicle, the demand for increasing helicopter usage for passenger 
transportation and the demand for high-speed maneuverability helicopters for defense 
have underlined the need for vibration reduction. Vibration reduction can be achieved in 
a number of ways. Vibration reduction schemes adopted in helicopters can be classified 
as either passwe or active control methodologies. The passive control scheme includes hub 
or blade-mounted pendulum alisorbers, anti resonant vibration isolation devices like dy- 
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namic anti resonant vibration isolator (DAVI), anti resonant isolation system (ARIS), and 
liquid inertia vibration eliminator (LIVE); structural modifications; and structural opti- 
mization. Active control methodologies include higher higher harmonic control (HHC), 
individual blade control (IBC), active flap control (AFC), and active control of structural 
response (ACSR). It may be noted that while HHC, IBC, and AFC control schemes are 
aimed at reducing the blade loads in the rotating frame, ACSR is employed in the non 
rotating frame to nullify the effect of vibratory hub loads on the fuselage. 

Control engineering is based on the foundation of feedback theory and linear system 
analysis. Accordingly, control engineering is not limited to any engineering discipline but 
is applicable to Aeronautical, Chemical, Mechanical, Environmental, Civil and Electrical 
Engineering. A control system is an interconnection of components forming a system con- 
figuration that will provide a desired system response. The basis for analysis of a system 
is the foundation provided by linear system, which assumes a cause-effect relationship for 
the components of a system. A component or process to be controlled can be represented 
by a block as shown in Fig 5.1. 

An open loop control system utilizes a controller or control actuator to obtain the 
desired response as shown in Fig 5.2. The open loop control system utilizes an actuating 
device to control the process directly without using sensors. 

A closed loop control system (Fig 5.3) utilizes an additional measure of the actual 
output to compare it with the desired output response. The measure of the output is 
called the fectlback signal. A feedback control system is a control system that tends to 
maintain a relationship of one system variable to another by comparing functions of these 
variables and using the difference as a means of control. As the system is becoming more 
complex, the interrelation of many controlled variables may be considered in the control 
scheme. 

Here, in this study we formulate a closed-loop (position feedback) control scheme (Fig 
5.4) for vibration minimization using the ACSR scheme. The concept of an ACSR scheme 
is based on the principle of superposition of two independent responses of a linear system, 
such that the total response is zero. In the case of helicopters, the fuselage is excited 
by the application of controlled external actuators at selected locations, such that the 
total response of the fuselage caused by the rotor loads and the external actuator forces 
is a minimum. The rotor loads are transmitted to the fuselage through the gearbox 
support structure. In this section, active vibration control of a highly idealised 2-D 
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system (Fig 5.5) of a helicopter is considered, Mp and Mu are the fuselage and hub 
masses respectively connected by a spring-damper-sensor-actuator system. The sensor 
and actuators are made of piezoelectric stack system, as shown in Figs. 5.6a and b. 

The actuator is attached between the hub and the fuselage through a mechanism as 
in Fig. 5.7. The relation between the horizontal and vertical deflection, force at the 
actuator to force transmitted to hub and fuselage structure is derived as follows. 


5.1 Formulation of relation between horizontal and ver- 
tical deflection 

The gearbox support mechanism before and after deformation is shown in Fig. 5.7 we 
have 


AO = \ a? - 


52 


( 5 . 1 ) 


From Fig. 5.7 it follows that 


AO = AO-f- 


(Vi - A'2) 




0,2 _ (5 + = (^2 _ ^ 


f H- AO(Xi - X2) 


52 bu>r. 


tip + AO{Xi - X2) 

4 4 


4 4 2 

Neglecting higher powers of u and (A"i — X 2 ) we get, 

huu- 


AO{X, - A2) = 


Hip 


But from Efi. 5.1 it follows that, 
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(5,2) 





Where, Uup corresponds to the axial tip deflection and (A'l - ^ 2 ) to the vertical 
deflection (Fig. 5.7). Thus the above relation 5.2 helps transform the deflection from 
sensor/actuator to that transmitted to the hub and fuselage. 


5.2 Formulation of relation between force at sensor/actuator 
and the force at hub/fuselage 

Let the force at sensor be tensile. The deformation is considered to be very small such 
that the angle T remains same for the deformed and undeforrned configuration. Then, 
the following cfinations are written considering each of the truss members to be axial 
bars. 


^Bc) ~ ( 5 - 3 ) 

(-^aZT "b ^Iw) ~ (5-4) 

where the subscripts (on LHS) denote the respective truss members and Fps denotes 
the force at sensor, whereas the force F<; denotes the force transmitted to hub/fuselage. 
Similarly, 


(^Ab ■+■ ^cd) cos 4/ - F PS 


(5.5) 


{Fm + Fpc) sin '^^Fs 


Considering Eqs. 5.3 and 5.5, 


^AB + 


DC 


Fad'^ Fnn 


CD 


From Eqs. 5.4 and 5.6, 


Fab "b F-, 


AD 


Fcd + F, 


BC 


(5.6) 


(5.7) 


(5.8) 
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Adding up Eqs. 5.7 and 5.8 we obtain, 

F-xb = ^CD = and similarly = F-jjy = 

Therefore, (Fi + Fj) cos T = Fps and (Fi + F 2 ) sin T = Fs 
i.e., tan^ = ^ 


Fs = Fps tan T (5.9) 

where from geometry tan Tis given as 

^ Jaq? - 1? 

tan = — = 

h 

The same formulation is carried out for actuation case, wherein Fps is substituted as 
FpA Fs is substituted as Fa- 

5.3 Analysis of Vibration Problem 

This study has been carried out in two sections, one in which only a sensor is provided 
Iretween hub and fuselage and the other includes both sensor and actuator. The former 
study enables us to find out the frequency response of the system. This frequency is then 
used as the excitation frequency in the formulation of latter part of the study wherein 
the actual vibration reduction problem using closed loop active vibration control scheme 
is solved. 

5.3.1 Sensing Mode 

The relation between sensing force, deflection, and potential induced (Fig. 5.6a) is given 
as (See Chap. 4) 


^tip — ttiF PS 


( 6 . 10 ) 


Z.C.j 


(5.11) 



and 


^induced — ^iP'pS (5.12) 

using Eqs. 5.2, 5.9, and 5.11 we get a relation between Fg and (X1-X2) in the form of 

Fs = -F\X, - X 2 ) (5.13) 

For given data as in Table 5.1 oi = 4.18031 * 10"^, cq = -7 = 2.392 * 10®, bi= 7.8955, 
and F* = 8.44 * 10® 

Formulation of equations of motion 

The equation of motion of the system in Fig. 5.5 can be written as 


MhXi + C'(A^ - .Y2) + K{Xr - X2) = Fs + Foe^^ 


Substituting for Fs, from Fq. 5.13 

MnX, + C{X, ~ X 2 ) + F(Xi - X 2 ) + F*(Xi - A' 2 ) = FoF^^^ 

This equation corresponds to hub/gear box motion. 

The equation corresponding to the fuselage motion can be written as, 


MpX2 + C(A2 - Xi) + K(X2 - Xi) + F*(yY2 - Xi) = 0 


Assuming a solution of the form Xi = A’lC'^'and A'’2 == A2e'“^ Then we have 

Ad 




iiLl 


re f 


^■o/A- /‘(ii)'' - (I + f )(i + i‘H±7r - mdtni + ^ 




and 


A, 


're f ' 


1 + ^ + i^ 

^re f 


‘^ref ^ 


An/A - (1 + f )(1 + - a((±;ni + /■) 


f ' 


^ref ' 


where. 


/.i = Mass Ratio 


(5.14) 


(5.15) 


Reference Frequency 
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Cref = 2MpOJref— Critical Damping Coefficient 


^ Damping Ratio 

The magnitudes of Eqs. 5.14 and 5.15 are plotted against the frequency ratio and the 
resonance/natural frequency can be identified from frequency response. 

5.3.2 Actuation mode 

Figure 5.6b represents an actuator. It is acted upon, by an external mechanical loading 
and voltage. The relation between the actuation force, applied voltage and deflection is 
given by equation (See Chap. 4), 

Uiip = CyFpA + CaRp/l 

FpA = CyUtip + C2VpA (5.16) 

_ VpA=G^induced (5.17) 

where G is the gain, ciand 62 arc constants. 

The gain is given as 

G = Ge"' 

Using Eqs. 5.2, 5.10, 5.12, 5.16, and 5.17 we get a relation between Fa and (X 1 -X 2 ) 
in the form of 


Fa = -F{Xy - X 2 ) 


(5.18) 


and when both sensor and actuator are included as two units in the system (Fig. 5.5) 
i.e., incorporating oq. 5.13 


F=-F{Xy-X 2 ) 
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where F is a linear function of gain and is given as 

F = di + daG 

where d.i = F + F* 

For given data as in Table 5.1 , ci = 5.6568G * 10^’’, C 2 = -1.87 * 10~®, F = 6.238 * 
10'’, ^2 = 2.2 * 10'’^ and hence F = 14.678 * 10''’ + 2.2 * lO^G. 

Formulation of equations of motion 
The equations of motion are given as 

MnXi + C{Xi - .Y2) + K{Xi - X2) = Fa + Fs + 

This corresponds to equation for hub/gear box, which can be written as 

MnX\ + C{Xi - X 2 ) + K(Xi - X 2 ) + F(Xi - X 2 ) = Foe^^ 
and the equation corresponding to the fuselage 


MfX 2 + C{X2 - Xy) + K{X2 - Xi) + F{X2 - X,) = 0 
Let A'] = A'le^^'and A :2 = Then we have 


+ F -. ( ,w. 1 


- (1 + f)(i + + 1^) 

V 1 + ^ + 7 ^ 

A 2 _ K ^ ^U),e/ 

- (1 + l)(i + ^‘)i±:r‘ - + /<) 


(5.20) 


The magnitudes of Eqs. 5.19 and 5.20 can be obtained for various values of gain G 
and phase rj. 

The response of a passive control system (Ref. [22]) are given as follows; 


Rh ^ .‘±^.2 \ 1 

Fz! K w + //uhaCl + L)^) — (1 + /.q + 112) — t 2 ^(^;^)(l + ;Ui + H2) 


88 



Rf f^ref~.2 


Pz/K 




OJ 


l-( 


+ L)) + zgg: 


“)^(/J'l + + /^l//'2(l + PV) ~ (1 + /^1 + (-‘‘2) — + Ml + M 2 ) 


where, 

Rn, Rp- Hub and Fuselage motions (Fig. 5.8) 

Pz -External force 

Mji, Mp, Mj - Hub, Fuselage and Isolator mass respectively 


Ml 

M 2 


Mr. 

M F 
Mj 
A4 F 

L = ^ 


5.4 Results and Discussions 

The frequency response of different systems is shown in Figs. 5.9. The Figs. 5.9 (a), 
(b), and (c) clearly depict the shift in the natural frequency of different systems (systems 
without sensor and actuator, effect of inclusion of sensor, and systems with both sensor 
and actuator). The reference frequency chosen was 6.78 Hz. It is observed that the 
natural frequency shifts from 19.5776 Hz (for a system without sensor and actuator), to 
21.5466 Hz for a system with sensor only and then to 22.8935 Hz for a system with both 
sensor and actuator. Figures 5.10 (a) and (b) show the phase and potential induced as a 
fuctioii of the frecpiency ratio, when the system is in sensing mode. Figures 5.11 depict 
the response of the system (in actuation mode) with respect to gain for different phase 
angles (of gain). It is observed from Figs. 5.11 (a)-(e) that a better vibration reduction 
can be obtained with increase in phase (from 0 to 180). Figure 5.11 (f) shows the response 
of the system when the damping of the system is increased to 0.04% (from the original 
value of 0.004%). It is observed that minimum fuselage vibratory response is observed at 
a gain value of 24.8 and phase of 180 deg. Figure 5.12 (a)-(c) show the frequency response 
of the actuator at various values of gain. Figure 5.13 shows the response of the system 
(Ref. [22]; Fig. 5.8) as a function of mass ratio. It is observed that minimum fuselage 
vibratory response is obtained when mass of the isolator is taken to be 1.3 kg. Using this 
value of mass of isolator the plots 5.14 are drawn. Figure 5.14 depicts the response of the 
system given in Ref. [22] when isolator mass is fixed at a value of 1.3 kg. Figures 5.15 
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and 5.16 give the response of tlie sensor and actuator when steel is chosen as the metallic 
layer. 

5.5 Summary 

A closed-loop active vibration control problem using position feedback is solved in this 
chapter. The results that were obtained from chapter 4 are used while writing the equa- 
tions of motion. Firstly the frequency response of the system is found out and then 
that utilized for solving the actual vibration problem. The results of this formulation 
were compared with the results given in Ref. [22 j and in latter case it was shown that the 
inclusion of an isolator betw^een the hub and the fuselage reduces the fuselage vibratory 
loads and increases the hub loads significantly, whereas our model simultaneously reduces 
both the fuselage and hub vibratory loads. This observation is of great importance in the 
industry. 
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Table 5.1: Data for vibration control 


Diameter of the bar = 2 rnm 

Length of the bar = 100 mm 

Number of layers of alternate MP...M = 101 

Thickness of metal layer = 0.5 rnm 

Thickness of piezo Ia.3fer = 1.5 mm 

Young’s modulus for steel = 208 GPa 

Young’s modulus for Al = 70 GPa 

M,,= 65 kg 

R = 6 m 

D= 32 rad/s 

K = M,/ 

, , — nn 

OJref — Y Mp 
Ci.ef=2*Mp*cu,e/ 

C=0.033 

Mass of fuselage = 2200 kg 
Mass of hub = 300 kg 
L=12.5 
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Figure 5.1: Process under Control 



Figure 5.2: Open Loop Control System (No Feedback) 


Actiicil 



Figure 5.3: Closed-Loop Feedback System 
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Figure 5.4: A Schematic Diagram of Position Feedback Control System 



S - Sensor 
A - Actuator 

Figure 5.5; Schematic of gear box support structure 
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(a) Sensor (b)Actuator 


Figure 5.6: Piezo stack sensor and actuator 


A 



(a) (b) 


Figure 5.7; Actuator model before and after deformation 
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Figure 5.8: Rotor luib/isolator/fuselage dynamical model (Ref. [22]) 
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(a)without sensor and actuator 




(b)with sensor only 




(c)with both sensor & actuator 


Figure 5.9: Variation of displacenieut and acceleration of different sytsems with frequency 
(metallic layer- Aluminium) 
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Phase of X 



(a) (b) 



Figure 5.10: Variation of phase and potential induced with frequency (sensing mode 
metallic layer- Aluminium) 










, '“SiofXK/Fo) . . log,o(XK/Fn) 




(a) Gain = -23 (i.e., 77=180^) 






Figure 5.12: Frequency response of the actuator at various values of gain (metalli layer 
Aluminium) 
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Chapter 6 

CONCLUDING REMARKS 


There has been a lot of research being undertaken to study the behavior of smart ac- 
tuators, their modeling and how they can be useful in real life applications. One such 
attempt is made in this thesis. We present electromechanical and thermal modeling of the 
piezo stack actuator. This includes modeling of hysteresis and butterfly loops, modeling 
of the actuator to obtain deflections and induced potentials, heat transfer analysis, and 
last but not the least is the analysis of the vibration problem. The important conclusions 
of this study are as follows:- 

1. The domain switching model for hysteresis and butterfly loops is a much preferred 
model as compared to the complex stiffness model as it generates both the loops 
even at small values of electric field. 

2. The inclusion of convection has a significant influence on the temperature profile of 
the piezo stack actuator. The temperature profile obtained in all cases is symmetric 
about the center of the bar. 

3. A linear relation is obtained between tip deflection and applied force in sensing; tip 
deflection, applied force, and applied voltage in actuation . It was observed that 
there exists a linear relation between potential or voltage applied and the block 
force (force applied to negate the tip deflection obtained through applying voltage). 

4. Using a piezo stack actuator as both sensor and actuator, a closed loop vibration 
control problem of a highly simplified 2-D model of a helicopter has been studied. 
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The results show that it is possible to reduce the vibratory response in both the 
hub and fuselage contrary to a passive vibration control system. 
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Chapter 7 

SCOPE FOR FUTURE WORK 


The present formulation can be considered as an initial step in the development of overall 
modeling of the smart actuators. It can lead a way to a more focused research in this 
area. A complete theoretical analysis has been presented over here. So as to validate 
this analysis an experimental setup can be established. Coming up to the vibration 
control part, this work can be extended to incorporate it in a real life helicopter vibration 
problem. 
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